Surface and bulk properties of ballistic deposition models with bond breaking 
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We introduce a new class of growth models, with a surface restructuring mechanism in which 
impinging "hot" atoms may dislodge suspended particles (deposited following the ballistic deposition 
rule), previously aggregated on the same column in the deposit. The flux of these "hot" particles 
is controlled through a probability p. These systems present a crossover, for small values of p, 
from random to correlated (KPZ) growth of surface roughness, which is studied through scaling 
arguments and Monte Carlo simulations on one- and two-dimensional substrates. We show that the 
characteristic time of crossover t x scales with p according to tx ~ P~ v with y = (n + 1) and that 
the interface width at saturation W sa t scales as W sa t ~ P~ S with S = (n + 1) /2, where n is either 
the maximal number of bonds breaked or of suspended particles that may be dislodged. Thus, the 
exponents y = 1 and 5 = 1/2 or y — 2 and 5 = 1 found in all previous works, are not universal. 
Using scaling arguments, we show that the bulk porosity P of the deposits scales as P ~ p y ~ s for 
small values of p. This general scaling relation explains previous results present in literature and is 
confirmed by our numerical simulations. 
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I. INTRODUCTION 

The morphology and evolution of the surface and 
porous structure in thin-film growth has been a subject 
of extensive theoretical and experimental interest [l], Q ■ 
From a theoretical side, a number of surface growth mod- 
els have beenproposed, being the ballistic deposition 
(BD) model Q one of the first. The BD model was 
initially formulated to explain sedimentary rock forma- 
tion, but has also been extensively studied as a model 
of low-temperature thin-film growth and surface rough- 
ening In the last decades, the surface properties of 
the BD model were widely studied (see refs. |l|, y, |j| and 
references therein) and it is believed to be in the Kardar- 
Parisi-Zhang (KPZ) Q universality class. More recently, 
some efforts were also devoted to understand its bulk 
structure [1, 0|- Furthermore, some variations of the BD 
model were also considered. Some examples are ballistic 
deposition with oblique incidence of particles @, Q , slip- 
pery BD 0, [HI EH , and competitive growth models in 
which particles are deposited follo wing pr obabilistically 
either BD or other deposition rules |12 - [20l |. 

In all models above, after deposition, the particles are 
permanently aggregated in the deposit and do not change 
their positions (limited mobility models) . Here, we study 
a different class of BD-like models in which aggregated 
particles can be rearranged in the deposit due to bond 
breaking, hereafter referred to bond-break ballistic depo- 
sition (BBBD) models. More specifically, BD particles 
are deposited with probability p and "random deposition- 
like" (RD-like) particles with probability 1 — p. When 
a "RD-like" particle impinges over suspended particles, 
their bonds can break and they fall together until reach- 
ing the next top of the column. We consider the case 
where up to n bonds can break (BBBDB models) and 



the case where up to n suspendend particles can change 
their positions, independently of the number of lateral 
bonds (BBBDS models). When bonds do not break, the 
"RD-like" particle is deposited on the top of the column 
(simple RD rule). 

In a very simplified way, a ballistic deposition mod- 
els a low-temperature system (with thermally activated 
diffusion process inoperative) where all impinging atoms 
immediately dissipate the energy released upon forma- 
tion of the atom-surface bond. However, in general, such 
energy could be not immediately dissipated and the im- 
pinging "hot" atom may perform/produces some tran- 
sient movements at surface. Some examples are down- 
ward funneling, transient mobility, knockout and cascad- 
ing knockdown of adatoms, which are processes respon- 
sible for smooth the surface in low-temperatures [2TI — 
I23I ]. In the last one, the "hot" atom dislodge previously 
deposited atoms that are only partially supported (sus- 
pended atoms), and they can move together to the "bot- 
tom valley" . Such mechanism induces the formation of 
a more compact film, eliminating overhangs at the sur- 
face. Thus, the bond break in the BBBD models may 
represents (in a simplified way) this concerted multiatom 
dynamics. 

Since for p = and p = 1 the BBBD models behave 
as plain RD and BD models, respectively, they present 
a competition from random to correlated (KPZ) growth. 
In systems with this kind of competition, the scaling be- 
havior of the roughness (or interface width) W is given 
by the empirical form [3, [2(| : 

where F(u) is a scaling function which behaves as F(u) ~ 
u Pbd with /3bd = zbd/o-bd for u << 1 and F(u) ~ 



const, for u » 1. olbd is the roughness, (3bd is the 
growth and zbd is the dynamical exponents of the BD 
model. In general, 8 and y are exponents dependent of 
the correlated model and independent of the dimensional- 
ity of t he sy stem. Results for several competitive models 
EH, Hfl 0] suggests that the exponents 8 and y are 
not directly linked to the universality class of the corre- 
lated process and can be separated in two groups: the 
first one including solid-on-solid (SOS) correlated mod- 
els, for which 8=1 and y — 2; and the second one includ- 
ing correlated models with lateral aggregation (ballistic- 
like models), for which 8=1/2 and y = 1. These pair of 
exponents are claimed to be universal, i. e., every model 
with a crossover from random to correlated growth would 
present one of those pair of exponents |20j. As we will 
show below, it is not the case, and for ballistic-like mod- 
els the exponents assume the general values 8 = (n+l)/2 
and y = n+ 1, where n is the maximal number of cither 
bonds breaked or suspended particles dislodged. 

Beyond the surface properties, in practical applications 
the internal structure of the material (beneath the sur- 
face) may be also relevant to determine mechanical and 
transport properties of the system. In this context, some 
recent works have focused also the porous (bulk) struc- 
ture [1, HH-Hll . Here, we study the bulk properties of the 
deposits generated by the BBBD models. Using scaling 
arguments, we show that the porosity P of the deposits 
scales as p v ~ s for small p, where 8 and y are the expo- 
nents defined above (Eq. [T|). Thus, this work advances 
over previous ones, which only relate P with p [ 1 lL 1 2 71 . 1 2 81 ] . 
since here we relate directly the scaling exponents of the 
amplitudes of surface roughness with that of bulk quan- 
tities. 

The rest of this work is organized as follows. In Section 
HI1 we describe the BD models with bond-breaking intro- 
duced here. In Section |TTT] we present the scaling theory 
to explain the behaviors of surface and bulk properties 
for small values of the probability p. Numerical results 
for the BBBDB and BBBDS are presented in Subsec- 
tions [IVA] and |IVBJ respectively. Final discussions and 
conclusions may be found in Section [V] 



II. THE MODELS 

In the BBBD models considered here, the impinging 
particle, following a vertical trajectory, can either dissi- 
pate all its energy at collision sticking at the point where 
it first contacts the surface (BD rule) or transfer energy 
to the particles bellow it, dislodging these particles if they 
are suspended. In this case, the particles (impinging and 
dislodged) are deposited at the next top bellow them. 
Two kinds of models were considered when dislodging 
particles: i) accounting the number of bonds breaking 
(BBBDB models); and ii) accounting the number of sus- 
pended particles dislodged (BBBDS models). In the first 
case, up to n bonds can break, and, in the second, up to 
n suspended particles can be dislodged. When it is not 
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FIG. 1: Illustration of the bond breaking rules (happen with 
probability 1 — p) for the BBBDB models with a) n — 1 and 
b) n = 2, and BBBDS models with c) n = 1 and d) n — 2. 
Deposits before (left) and after (rigth) deposition of seven 
particles are shown. 



possible to dislodge particles, the impinging particle is 
deposited at the top of the incidence column (RD rule). 

In both classes of models, a particle is dropped ver- 
tically at a randomly selected column i and sticks on 
the first site encountered on the surface that is nearest- 
neighbor of an already deposited particle (BD rule) with 
probability p. In BBBDB models, with probability 1 — p, 
the number k of bonds of the target (top most) particle 
at column i is determined. If k > n the incident parti- 
cle is deposited at the top of the column, without bond 
break. If k = n the target particle is unbounded and slide 
down vertically (together with the incident particle) to 
the next local minimum. Note that, in most cases there 
is no movement of the target particle because it is not 
suspended (there is a particle just bellow it). If k < n, 
particles bellow the target one could be unbounded (and 
dislodged also), provided that the total number of bonds 
breaking are smaller or equal to n. Figs. Q] (a)-(b) illus- 
trate this growth rule for n = 1 and n = 2. 

In BBBDS models, with probability 1 — p, the num- 
ber k of consecutive suspended particles in column i is 
counted, starting from the top. If k > n, the incident 
particle is deposited at the top of the column, without 
bond break. If k = n the block of n suspended particles 
and the incident one slide down vertically to the next top. 
If k < n, suspended particles bellow the fist suspended 
block of particles may also be dislodged, until a block 
of size smaller or equal n be attained. Figs. [1] (c)-(d) 
illustrate this growth rule for n = 1 and 2. 
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For n = 0, BBBDB and BBBDS models reduce to the 
simpler RD-BD model proposed by Horowitz and Albano 
flij . Thus, the RD-BD model is a particular case of the 
BBBD models, referred here just as BBBDO model. Sim- 
ulations of this model showed that it exhibits a crossover, 
at a characteristic time t x , from RD to KPZ growth and 
gave 8 ~ 1/2 and y ~ 1 in d = 1 + 1, d = 2+1, and 
d = 3 + 1 dimensions (see also [2(| and references 
therein). For completeness, we present bellow simula- 
tions results of the BBBDO model, that are consistent 
with the results above. Moreover, we present also re- 
sults for the bulk properties of this model, which was not 
studied yet. 



III. SCALING THEORY 

Different scaling arguments have already explained the 
values of the exponents 8 and y obtained in simulations 
of several competitive lattice models jTR 21) . For the 
BBBDO (RD-BD) model, the basic argument is that, for 
small values of p, the average time for a correlated BD 
event (which involves lateral aggregation with probability 
p) to take place at a given column is r ~ 1 /p. During the 
period r, particles on average are directly deposited onto 
the surface according to the simple RD rule and the local 
height increases by y/r. Thus, the BBBDO model can be 
viewed (for small p) as a limiting BD model in which 
time scales as r ~ 1/p while roughness scales as W ~ 
\fr 1/p 1 / 2 . This explains the scaling function given by 
Eq. ([1]) and the conjectured values 8 = 1/2 and y = 1 
for this model. This argument also gives the relationship 
S = y/2 which is valid for all models with crossover from 
RD to some correlated deposition dynamics [2Qj . 

For others BBBD models, we expect a similar argu- 
ment. However, due to the bond breaking processes, 
which is frequent for small p, a single BD event is not 
sufficient to introduce lateral aggregation and cancel the 
random fluctuation of the heights. In order to ensure 
lateral growth, it is necessary that more BD particles 
fall (consecutively) at the same column. For example, if 
one suspended particle could be dislodged in the deposit, 
two consecutive BD events are needed to produce a cor- 
related growth, giving r ~ 1/p 2 and W ~ y/r ~ 1/p, 
and the exponents 8=1 and y = 2. In general, if n 
bonds could break, or n suspended particles could be 
dislodged, it is necessary n + 1 ballistic particles to in- 
troduce correlations, which leads to t ~ l/p™ +1 and 
W ~ ^fr i/p("+ 1 )/ 2 j and thus to the expected ex- 
ponents 



5 = (n + l)/2 



and 



y = {n + 1) 



(2) 



for these models, in any substrate dimension. It is worth 
mentioning that in all previously studied competitive 
models involving RD and solid-on-solid (SOS) correlated 
models (for example, RD-RSOS (restricted SOS) [19J and 
RD-RDSR (random deposition with surface relaxation) 
[20j|). exponents 8 — 1 and y = 2 were found. In contrast, 



systems with RD and BD-like models (for example, the 
BBBDO model [15|,|20j]) have the exponents 5=1/2 and 
y = 1. This resulted in a classification of the models in 
two groups: RD-SOS models (with 8=1 and y = 2) 
and RD-BD like models (with 8=1/2 and y = 1) 
[lil [l9l I20I [24]], and these two sets of exponents were 
claimed to be universal (the only two pair of possible 
values) in random-to-correlated crossover (20j . Our re- 
sults shows that for RD-BD like models, the exponents 
would be 8 = (n+ l)/2 and y = n + 1, that gives 8=1/2 
and y = 1 in the case n = (no bond breaking) already 
considered. 

It is possible to extend the scaling arguments presented 
above to explain the behavior of the bulk porosity P of 
the deposits, for small p. The porosity P is the fraction 
of empty lattice sites inside the deposit: 



P = 



V P 



Vs + Vp 



(3) 



where Vp is the total volume of the pores (vacant sites) 
and Vs is the total volume of the solid (particles) in the 
deposit. For small p, the pores are narrow in the horizon- 
tal direction (~ 1 site width) but vertically high. A typ- 
ical height difference between neighboring columns is of 
order yfr ~ p~ s for small p and, thus, the lateral aggrega- 
tion events create pores with volume of order p . After 
Vs depositions, the average number of pores created is of 
order p y Vs, since the characteristic time in which corre- 
lations between neighboring columns are built scales as 
t ~ p~ v . Thus, we obtain P ~ p A ? with A = y — 8, 
which must be valid in any substrate dimension. It is 
worth mentioning that for the bidisperse BD model (for 
which 8=1/2 and y = 1) the scaling behavior P ~ p 1 / 2 
was obtained in Ref. [27|, consistent with the most gen- 
eral scaling above. Here, we must expect the exponents 
A = (n + l)/2 for the BBBD models. 



IV. NUMERICAL RESULTS 

The following results are for simulations on d = 1 + 1 
substrates of sizes in the range L = 64 — 1024, and d = 
2+1, with sizes in the range L = 8 — 128. In both cases, 
we assume periodic boundary conditions. Averages over 
400 — 1000 different runs were made. 



A. BBBDB models 

Figure [2] shows plots of the roughness W versus time 
t obtained for the BBBDB 1 model for two values of p, 
in d = 1 + 1. For short times, say t < t x , the growth is 
dominated by the RD process and W ~ t 1 / 2 for p < 1. 
At intermediate times, t x < t < ts, correlations have 
developed and the BD process dominate giving W ~ t?. 
Using data for different lattice sizes and extrapolating to 
L — > 00, we obtained growth exponents /3 ~ 0.33 for ev- 
ery studied value ofp>0inrf=l + l. This procedure 
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FIG. 2: Roughness W versus time t for the BBBDB1 model 
with p = 0.05 and p — 0.5, for a system size L — 512. The 
inset shows effective growth exponents /?l against (l/L) ' 25 
for p = 0.5. The linear fit gives the extrapolated value — 
0.35 ±0.03. 

to obtain /?l^oo is illustrated in the inset of Fig. [21 for 
p = 0.5. This asymptotic exponent is consistent with the 
expected KPZ one (f3 = 1/3) and with previous simula- 
tion results for the BD model in d = 1 + 1 [4] . Finally, for 
t > tg the correlations can no longer develop due to the 
geometrical constraint of the lattice size and saturation 
occurs. Similar W x t behaviors are found in all models 
studied here and in other systems with random to KPZ 
crossover [2fil [27j. 

Bellow, we present results for the BBBD0 and BBBDB 
models with n — 1 and 2, in d = 1 + 1 and d = 2 + 1 
substrates. 

In Fig. [2a) the saturation roughness W sa t as a func- 
tion of p is shown. As expected, for a given p, a larger 
n implies a larger W sa t- For small values of p, straight 
lines are observed in the log- log plot, in agreement with 
Eq. fl} , and the best fits give the exponent 6 in excellent 
agreement with the scaling theory of Sec. IIIII (see Table 

The crossover times t x , for different values of p, are 
shown in Fig. |3fb). The straight lines in the log- log 
plot, for small values of p, show that the scaling behavior 
tx ~ P y from Eq. [T]is satisfied and yields exponents y in 
great accordance with the predicted ones, as summarized 
in TableU 

As the roughness, the porosity of the aggregate grows 
in time and attain a stationary value P for large t. Fig. 
[3Uc) shows this stationary porosity as a function of the 
probability p. A good scaling behavior is found for small 
p and the exponents A obtained are consistent with the 
expected values (see Table [J). 

The results above are for substrate sizes L = 128 (in 
d = 1 + 1) and L = 32 (in d = 2 + 1). We notice 
that we have to work with small system sizes because the 
saturation times for p <^ 1 became very long for large n. 
However, it is not a limitation since, as pointed out in 
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FIG. 3: a) Saturation roughness W 3a t, b) crossover times t x 
and c) porosity P versus probability p for the BBBDB models 
with n = (red triangles), n = 1 (blue squares) and n — 2 
(green diamonds). Data for simulations on d — 1 + 1 sub- 
strates of size L — 128 (open symbols) and d — 2 + 1 with 
L — 32 (full symbols) are shown. In order to improve visual- 
ization, in (a) and (c) the data for n — 1 and 2 in d = 2 + 1 
are, respectively, multiplied (shifted up) and divided (shifted 
down) by 1.5. The straight lines have slopes (from the bottom 
to top): a) —1/2, —1 and —3/2, b) —1, —2 and —3, and c) 
3/2, 1 and 1/2. 



Ref. [20|| , the exponents <5 and y have weak finite-size 
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BBBDO 


BBBDB1 


BBBDB2 


y 


0.98(3) 


1.96(5) 


3.0(1) 


6 


0.46(5) 


0.98(3) 


1.50(2) 


A 


0.43(3) 


0.97(4) 


1.48(4) 


y 


1.02(5) 


2.02(3) 


3.03(7) 


6 


0.47(4) 


0.98(3) 


1.49(3) 


A 


0.40(5) 


0.96(5) 


1.50(3) 



TABLE I: Amplitude exponents y and S and porosity expo- 
nent A = y — S for d = 1 + 1 (top) and d = 2 + 1 (bottom) . 



effects. In fact, from simulations for other system sizes 
(smaller and larger in same cases) we obtain exponents 
close to that shown in Table fl] 

Summing up, our numerical results are consistent with 
the roughness scaling relation (JXJ) , and the exponents 5, 
y and A obtained confirms the exact values from the 
scaling theory presented above. 



B. BBBDS models 

We made simulations of BBBDS models with n = 1, 
2 and 3 in d = 1 + 1 and d = 2 + 1 dimensions. As 
discussed above, the roughness evolution with time has 
a similar behavior to the one shown in Fig. [2] for all 
BBBDS models with p < 1. 

In BBBDS models, more particles are dislodged, for a 
given n, than in BBBDB ones (see Fig. [TJ. Thus, in the 
former the correlations develop more slowly, what makes 
the crossover times and the saturation roughness larger. 
On the other hand, the porosity is smaller in BBBDS 
models, since more particles are compactified by dislodg- 
ment. For example, for n = 1 and fixed L and p, we found 

t BBBDS /t BBBDB ^ Lg m d = 

1.2 



BBBDB ^BBBDij j^Bi 

1 + 1 and ps 3 in d = 2+1, while P bbbdb / p bbbds 
in d = 1 + 1 and ps 1.5 in d = 2 + 1. 

In Fig. @]we show the saturation roughness W sa t, the 
crossover time t x and the stationary porosity P as func- 
tions of p. As in BBBDB models, we found here good 
scaling behaviors, for small p, consistent with the scaling 
relation given in Eq. [TJ The exponents y, S and A are 
in great agreement with the expected ones, as shown in 
Table HI 

The data above are for substrate sizes L = 128 in 
d = 1 + 1. In d = 2 + 1 we have L = 32 for n = 1 
and 2, and L = 16 for n — 3. Exponents obtained for 
other system sizes are close to the ones in Tab. UH again, 
showing that such exponents have negligible finite-sizes 
effects. For n = 3, the large crossover and saturation 
times restrains our simulations to substrate sizes L < 16 
in d = 2 + 1. For these sizes, there is not a clear KPZ 
region in the W xt curves. In other words, the crossover 
and saturation times became very close, so that it is not 
possible to determine the former one. Thus, we are not 
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FIG. 4: a) Saturation roughness W 3a t, b) crossover times t x 
and c) porosity P versus probability p for the BBBDS models 
with n = 1 (red triangles), n = 2 (blue squares) and n = 3 
(green diamonds). Data for simulations on d — 1 + 1 sub- 
strates (open symbols) of size L — 128 and d = 2 + 1 (full 
symbols) with L — 32 (for n = 1 and 2) and L — 16 (n = 3) 
are shown. In order to improve visualization, in (c) the data 
for d = 2 + 1 are divided (shifted down) by 1.4. The straight 
lines have slopes (from the bottom to top): a) — 1, —3/2 and 
-2, b) -2, -3 and -4, and c) 2, 3/2 and 1. 



able to obtain the exponent y directly from the t x scaling. 
However, since A = y — 5, from the A and S exponents 



6 





BBBDS1 


BBBDS2 


BBBDS3 


V 


2.01(4) 


3.04(5) 


4.0(1) 


8 


0.99(1) 


1.50(2) 


2.0(1) 


A 


0.98(4) 


1.49(2) 


1.98(3) 


y 


2.1(1) 


3.05(5) 




5 


0.97(4) 


1.50(2) 


1.93(5) 


A 


0.99(3) 


1.52(2) 


2.02(4) 



TABLE II: Amplitude exponents y and 5 and porosity expo- 
nent A = y — S for d = 1 + 1 (top) and d — 2 + 1 (bottom). 



we obtain y s» 3.95. 

V. DISCUSSIONS AND CONCLUSIONS 

We have studied ballistic-like models where bond break 
of suspended particles is allowed (the BBBD models). 
The competition between RD-like and DB-like deposi- 
tions in these models produces a crossover from random 
to KPZ growth in the temporal evolution of the surface 
roughness that follows the same scaling relation observed 
in previous models with random to correlated growth. 
Using scaling arguments we show that the exponents de- 
scribing the roughness and crossover time amplitudes, for 
the BBBD models, are 8 = (n + l)/2 and y = (n + 1), 
where n is either the maximal number of bonds breaked 
or of suspended particles dislodged. Such exponents are 
confirmed by numerical simulations of the models on one- 
and two-dimensional substrates (sumarized in Tables U 
and[Tl| and they are expected to hold in any dimension. 
These results generalize the classification present in liter- 
ature, where models with random-to-correlated competi- 
tion were divided in two groups: solid-on-solid ones, with 



S = 1 and y = 2, and ballistic-like ones, where 8=1/2 
and y = 1. In fact, this exponents are not universal, as 
claimed in previous works, and in ballistic-like models 
we must have the general exponents found here (Eq. [2]). 
These exponents reduces to S = 1/2 and y = 1 if n = 0, i. 
e., if there is no bond breaking, as considered in previous 
works. 

Models accounting for the maximal number of bonds 
breaking (BBBDB) and for the maximal number of sus- 
pended particles dislodged (BBBDS) have different rates 
of particle dislodgement, what leads to large differences 
in crossover times, saturation roughness and porosity. 
However, they have the same exponents (given by Eq. 
[5]), showing that these ones are independent of some de- 
tails of the model and must be observed when n particles 
can be dislodged in a column. In fact, in Fig. [T]we see 
that in both classes of models up to n suspended parti- 
cles in a column could be dislodged, but in the BBBDB 
models these particles sometimes do not move because 
they have lateral neighbors, and thus a large number of 
bonds. 

Finally, we showed that the bulk porosity P of the 
deposit is related with the probability p (for small p) 
through the surface exponents discussed above, with P ~ 
py- s _ This scaling relation was confirmed with numerical 
simulations on d = 1 + 1 and d = 2 + 1 substrates, and 
must be valid in any substrate dimension. It also explains 
the relation P ~ p 1 / 2 found for the bidisperse BD model 
in [27j . and it must works in any system with random- 
to-correlated crossover forming a porous deposit. 
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